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Evolution of Primordial Magnetic Fields from Phase Transitions 
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We consider the evolution of primordial magnetic fields generated during cosmological, electroweak 
or QCD, phase transitions. We assume that the magnetic field generation can be described as an 
injection of magnetic energy to cosmological plasma at a given scale determined by the moment 
of magnetic field generation. A high Reynolds number ensures strong coupling between magnetic 
field and fluid motions. The subsequent evolution of the magnetic field is governed by decaying 
hydromagnetic turbulence. Both our numerical simulations and a phenomenological description 
allow us to recover "universal" laws for the decay of magnetic energy and the growth of magnetic 
correlation length in the turbulent (low viscosity) regime. In particular, we show that during the 
radiation dominated epoch, energy and correlation length of non-helical magnetic fields scale as 
conformal time to the powers —1/2 and +1/2, respectively. For helical magnetic fields, energy and 
correlation length scale as conformal time to the powers —1/3 and +2/3, respectively. The universal 
decay law of the magnetic field implies that the strength of magnetic field generated during the QCD 
phase transition could reach ~ 10~ 9 G with the present day correlation length ~ 50 kpc. The fields 
generated at the electroweak phase transition could be as strong as ~ 10" 10 G with correlation 
lengths reaching ~ 0.3 kpc. These values of the magnetic fields are consistent with the lower bounds 
of the extragalactic magnetic fields. 

PACS numbers: 98.70.Vc, 98.80.-k 
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I. INTRODUCTION 



Astronomical observations show that galaxies have 
magnetic fields with a component that is coherent over 
a large fraction of the galaxy with field strengths of the 
order of 10~ 6 G; see Refs. [l|-y| and references therein. 
Understanding the origin of these fields is one of the 
challenging questions of modern astrophysics. Generally 
speaking, there are two popular scenarios. The first one 
envisages the generation of magnetic fields through dif- 
ferent astrophysical mechanisms. More precisely, it is 
assumed that an initially tiny magnetic field is produced 
through a battery mechanism [|| ■ The correlation length 
of such a field is limited by galactic length scales. The 
second scenario to explain the origin of the magnetic field 
in galaxies and clusters presumes that the observed mag- 
netic fields were amplified from cosmological weak seed 
magnetic fields [5| . In this case the correlation length of 
such a seed field might be as large as the horizon scale to- 
day if we admit that the field has been generated during 
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inflation. There are different possibilities for seed mag- 
netic field amplification ranging from a magnctohydro- 
dynamic (MHD) dynamo to the adiabatic compression 
of magnetic field lines during structure formation 0, 0] ■ 

Galactic magnetic fields are usually measured through 
the induced Faraday rotation effect @, E| and, as men- 
tioned above, the field magnitude is of the order of a few 
10~ 6 G at typical scales of 10 kpc. The primordial mag- 
netic energy density contributes to the radiation field, 
and thus the big bang nucleosynthesis (BBN) bound im- 
plies fifi/ip < 2.4 x 10~ 6 [1] if the magnetic field has been 
generated prior to BBN. If the correlation length of the 
magnetic field is much larger than A# > 1 Mpc, smaller 
limits on the magnetic field energy density arise from the 
cosmological data, making i? max < a few 10 -9 G; see @ 
and references therein. A correlation length-dependent 
lower limit on magnetic fields in the intcrgalactic medium 
(IGM) could be derived from gamma-ray observations of 
blazars In the limit of large correlation lengths, 

\b > 1 Mpc, the bound is at the level of 10 _17 G (sec 
also Ref. [14[ for a discussions on possible uncertainties 
in the measurements of blazar spectra). 

It is possible, in principle, that the IGM magnetic 
fields are of primordial origin. Another possibility is that 
the fields are spread through the IGM by outflows from 
galaxies at late stages of the evolution of the Universe. 
To distinguish between these two possibilities, it is im- 
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portant to identify measurable characteristics of the IGM 
magnetic fields which are different in the two cases and 
to study the possibility of measuring such characteristics. 

In this paper we consider the observational properties 
of IGM magnetic fields expected if the fields originate 
from cosmological phase transitions (PT) such as the 
electroweak (EW) and QCD PTs (HG3. We follow the 
evolution of fields from the epoch right after the magne- 
togenesis up to the present day epoch. Our approach is 
different from that adopted in the previous studies, which 
mostly concentrated on the analysis of the range of pos- 
sible field strengths at a pre-defined scale of interest (e.g. 
IMpc). Instead, we study the evolution of field char- 
acteristics that are most relevant for measurements us- 
ing radio and gamma-ray astronomy. Specifically, we are 
interested in the evolution of the magnetic energy den- 
sity pb, which determines the characteristic field strength 
£?( cff ) = y/SnpB, and the characteristic correlation scale 
(integral scale) Am at which the field reaches the strength 
#( eff ). 

The present day integral scale depends on the tempera- 
ture T* and the number of relativistic degrees of freedom 
g+ at the moment when the primordial magnetic field is 
generated. These parameters determine the maximal al- 
lowed value of the magnetic energy density injected in the 
PT plasma, as well as the initial correlation length of the 
magnetic field [lj}. We do not separately consider the 
effect of helicity transfer related to the chiral anomaly, 
which might be important in the presence of strong mag- 
netic fields at temperatures above 10-100 McV (l8|. Such 
a transfer could be considered as part of the magnetoge- 
nesis process which could persist all the way down to 
the temperature scale of the QCD phase transition. We 
only use fundamental physical laws, such as conservation 
of energy, and how the magnetic field interacts with the 
cosmological plasma through MHD turbulence, and do 
not make any assumption about the physical processes 
responsible for the primordial magnetic field generation. 

In Sec. II we give an overview of the spatial and tempo- 
ral characteristics of the primordial magnetic field. The 
results of our analysis are presented in Sec. Ill, where 
we discuss the evolution of the magnetic field. Conclu- 
sions are presented in Sec. IV. We employ natural units 
with h = 1 = c and gaussian units for electromagnetic 
quantities. 



II. MODEL DESCRIPTION 



A. Effective Magnetic Field Characteristics 



We assume that the phase transition-generated mag- 
netic fields satisfy the causality condition [l6|, [2(| |21| . 
The maximal correlation length £ max for a causally gen- 
erated primordial magnetic field cannot exceed the Hub- 
ble radius at the time of generation, H~ . Hence 7 = 
timnx/H^ 1 < 1, where 7 can be associated with the num- 
ber of primordial magnetic field bubbles within the Hub- 



ble radius, N cx 7 s . The comoving length (measured 
today) corresponding to the Hubble radius at generation 
is inversely proportional to the corresponding PT tem- 
perature 
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(1) 



and is equal to 0.5 pc for the QCDPT (with g+ = 15 and 
= 0.15 GeV) and 6 x 10~ 4 pc for the EWPT (with 
<?* = 100 and T+ = 100 GeV), and the comoving primor- 
dial magnetic field correlation length £ max < Ah- This 
inequality assumes only the expansion of the Universe 
without accounting for MHD turbulence (free turbulence 
decay or an inverse cascade if a helical primordial mag- 
netic field is present). 

The maximal value of the primordial magnetic field en- 
ergy density must satisfy the BBN bound, i.e. the total 
energy density of the primordial magnetic field at nu- 
cleosynthesis ps (where a should not exceed 10% of the 
radiation energy density /9 ra d(aN)- In any case, the mag- 
netic field is generated by mechanical motion of charged 
particles, so that its energy density could constitute only 
a fraction of matter energy density, which is at most com- 
parable to the radiation energy density in the radiation 
dominated Universe. Note that the maximal value of the 
effective magnetic field is independent of the temperature 
at generation T*, and depends only very weakly on the 
number of relativistic degrees of freedom at the moment 
of generation. 

In what follows we are mostly interested in the evo- 
lution of the energy density of the magnetic field and 
the length scale which gives the dominant contribution 
to the energy density (the "integral scale") during the 
course of cosmological evolution. Taking this into ac- 
count we adopt the following idealizing approximation. 
We generate an initial primordial magnetic field by solv- 
ing the MHD equations for a certain time during which 
an external electromagnetic force is applied that is pro- 
portional to a delta function that peaks at the character- 
istic scale ko = 2tt/^ 1 . This corresponds to a magnetic 
field with correlation length £o- In this approximation 
the characteristic magnetic field strength at the scale £0 
is _B( cff ) = y/8irpB- We justify our assumption that £0 
should be identified with the size of the largest magnetic 
eddies by noting that the primordial magnetic field is 
involved in MHD processes driven by turbulence. It is 
natural to assume that the typical length scale of the 
magnetic field generated during the PTs is determined 
by the PT bubble size. 

The dynamical evolution of the coupled magnetic field- 
matter system leads to a spread of magnetic field over 
a range of scales. The resulting power spectrum of 
the magnetic field at small wavenumbers (or, equiva- 
lently, large distance scales) has the form of a power-law 
P M (k) = E A ,{{k)/(ATTk 2 ) = Ak UB with a normalization 



3 



constant A and a slope 1 ub- In particular j a white noise 
power spectrum corresponds to n# = [201, while the 
Batchelor spectrum corresponds to ns = 2 [2l[. The 
power law spectrum extends up to a time-dependent in- 
tegral scale above which the power contained in the 
magnetic field decreases rapidly due to turbulent decay 
and/or viscosity damping. 

Several previous studies, see [HI and references 
therein, describe the primordial magnetic field in terms 
of a smoothed (over length scale £) magnetic field B^ 
with B'l = (Bi(x)Bi(x)}\ ( . Knowing B^ and the slope 
of the power spectrum, one can calculate the strength 
of the smoothed magnetic field at any scale of interest £ 



B t = 



j3 (cff Vr(n B /2 + 5/2) 
(eAA/) ( " B+3)/2 ' 



(3) 



The smoothed magnetic field might be of interest in the 
context of certain problems. For example, the strength 
of a magnetic field smoothed over a scale £ ~ 1 Mpc is 
considered to be relevant in the context of seed magnetic 
fields for galactic dynamos. It is, however, important to 
note that B^ is strongly ns dependent for a given value 
of £. In particular, for causally-generated magnetic fields 
with ub > there is significant "magnetic power" only 
at small scales, and for £ ~ 1 Mpc the value of £?im P c 
is extremely small (25|. At the same time, it does not 
imply that the magnetic field itself is weak. In fact, it 
could be as strong as i?( eff ' ~ 10 -6 — 10 -7 G, close to the 
bound imposed by the BBN [I?} ■ Only in the case of a 
scale invariant magnetic field with ub — > —3, generated 
for example during inflation j26j , B^ is independent of £ 
and n_e, and is equal to B^ eS \ 



B. Phenomenological description of the magnetic 
field decay in the free turbulence regime 

After generation, the evolution of the primordial mag- 
netic field is a complex process affected by MHD as well 
as by the expansion of the Universe [27H331 [H, [3(| . In 
our description, to account for the expansion of the Uni- 
verse we make use of the fact that conformal invariancc 
allows for a description of MHD processes in the early 
Universe by simply rescaling all physical quantities in 
terms of their comoving values and using the conformal 



time r\ [T(|. After this procedure the MHD equations in- 
clude the effects of the expansion while retaining their 
conventional flat spacetime form. 

The magnetic evolution process strongly depends on 
initial conditions, as well as on the physical conditions of 
the primordial plasma. We need to determine the scal- 
ing laws for the following magnetic field characteristics: 
(i) magnetic energy density, (ii) correlation length, and 
(iii) magnetic hclicity. The magnetic energy and mag- 
netic helicity spectra are related through the realizabil- 
ity condition, \H M (k,7j)\ < 2E AI {k,7j)/k [H]. For the to- 
tal magnetic energy = J Enj(k,rj) dk and helicity 
UmW) = J H M (k,rj)dk we get 



where 



i M {ri) = 8^{ri) I Euik^k^dk 



(4) 



(5) 



is the comoving magnetic eddy correlation length (which 
corresponds to the physical integral scale Xm), initially 
set by the temperature at the magnetic field generation 
moment ^M,in = Ao = 7A# t ; see Eq. (TTJ), and is indepen- 
dent of the presence of magnetic helicity. 

Both helical and non-helical magnetic fields experience 
large-scale MHD decay resulting in an increase of the 
correlation length with a corresponding decrease in the 
ma gnet ic energy density at large scales; for a review see 
p9Tl37| . The time rate of this process depends strongly 
on the presence of magnetic helicity [13, . Taking this 
into account, we consider the cases of helical and non- 
helical magnetic fields separately. 

As noted above, the initial magnetic field configuration 
is given by a sharply peaked spectral energy density. The 
coupling between primordial magnetic field and plasma, 
which ensures the spreading of the fixed scale primor- 
dial magnetic field over a wide range of length scales, 
forms a modified magnetic field spectrum within a few 
turnover times (see Ref. [23[ for more details). The fi- 
nal realization of the spectrum is given by the Batch- 
elor spectrum, ns = 2. Our numerical simulation re- 
sults are in perfect agreement with the "causality con- 
straint" that in the cosmological context has been dis- 
cussed in Ref. 16 1, and studied through analytical ap- 
proach in Ref. [2l| 2 so that the power contained in the 



large-scale modes is very small, while the total magnetic 
energy density is sufficiently large. MHD processes also 



1 In general the magnetic field spectrum is determined through the 
Fourier transform (k) of the two-point correlation function of 
the magnetic field, (Bi(x)Bj(x + r)), with the spectral function 

El 



F^ I (k) = P ij (k) 



E M (k) 
Airk 2 



H M {k) 
Snk 2 



(2) 



Here Pij(k) 



kikj/k 2 , Eiji is the antisymmetric tensor, 



E (k) and H (k) are the magnetic energy and helicity spectra. 



2 A recent study based on semi-analytical calculations |22| showed 
the same shape for phase transition-generated magnetic fields. 
In laboratory plasma as well as in numerical simulations, the 
resulting magnetic field spectrum at large scales can be given by a 
white noise spectrum, ub = or even by a flatter Kazantsev 
spectrum, ng = — 1/2 |29|| . However, here we consider the case 
of a cosmological magnetic field for which the correlation length 
is strongly limited by the Hubble horizon. 
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are responsible for the generation of fluid perturbations 
when an initial magnetic field is present. These processes 
finally result in equipartition between magnetic and ki- 
netic energy densities [HI, [H[ . In contrast to the mag- 
netic field, the velocity field has a white noise spectrum, 
i.e. P K = E K /(4nk 2 ) = A K k nK with n K = due to 
the possible presence of longitudinal modes. To describe 
adequately the evolution of fluid motions coupled to the 
magnetic field we need to solve the complete set of MHD 
equations; see below. 

Below we present a phenomenological description of 
the MHD decay laws at large scales. In the case of he- 
lical fields we mostly follow the description presented in 
Sees. 4.2.3 and 5.3.2 of Ref. [H and Sec. 7.3.4 of [39|. 
On scales below £mj magnetic power is transferred to 
smaller scales via the so-called direct cascade by turbu- 
lence until it is finally damped at the smallest scale 
In MHD the magnetic field damping is usually deter- 
mined through the Reynolds number as £,d/£,M = i?e~ 3 / 4 
[37l |. The kinetic and magnetic Reynolds numbers in 
the early Universe can be extremely high, and thus one 
may expect that ^ <C £m- In both helical and non- 
helical cases the dissipative region of the energy den- 
sity spectrum is given by a Kolmogorov-type spectrum 3 



E 



M 



r< 2/3,. 

C K £m k 



5/3 



where Ck is a constant of order 



unity (1.6-1.7 for a wide range of Reynolds numbers), 
Em is the magnetic energy dissipation rate per unit mass, 
given by em = 2A k 2 E]\,,i(k) with A being the mag- 
netic diffusivity. At large scales above £m the magnetic 
field decay is strongly dependent on the presence of mag- 
netic helicity. The high conductivity of plasma ensures 
magnetic helicity conservation that is responsible for the 
transfer of spectral energy from small to large scales via 
a so-called inverse cascade. In the non-helical case the 
process is more complicated. As we will sec below, mag- 
netic helicity conservations leads to a faster growth of the 
correlation length and a slower decay of total magnetic 
energy. 



1. Non-helical magnetic fields 

As we already stated above, causality requires that 
EM{k) oc k 4 and this is a consequence of the divergence- 
free condition for the magnetic field [21| . On the other 
hand, there is no zero-divergence requirement for the ve- 
locity field, and this allows for the possibility to have a 
white noise spectrum for the velocity field, i.e. Ex{k) oc 
k 2 : see [2(|. We would like to note that our numerical 



simulations allowing the longitudinal forcing, see j35ll38l|. 
show a white noise spectrum for the velocity field as a fi- 
nal configuration. Under these conditions, the power of 
magnetic field modes on the large scales is much smaller 
than the power of plasma motions. Thus, potentially the 
magnetic field might be amplified via a transfer of energy 
from plasma motions at large scales. The time scale r\ on 
which the field can be amplified at large scales can be 
deduced from the induction equation 



dB v K {L)B 
~dt L 



77 ~ L/vl 



(0) 



i.e., the characteristic field amplification time scale is ap- 
proximately the plasma eddy turnover time. On this time 
scale an equipartition between kinetic and magnetic en- 
ergies could be reached over the distance range L. In the 
final configuration the equipartition between magnetic 
and kinetic energies is a consequence of the coupling be- 
tween magnetic and velocity fields. 4 The growth of the 
magnetic field up to equipartition with the fluid on large 
scales is somewhat similar to the phenomenon of an "in- 
verse cascade" of the magnetic power spectrum. Note 
however that the source of power in the inverse transfer 
of non-helical magnetic fields is different from the power 
source in the case of non- helical fields. In the case of non- 
helical fields, the power in the large wavelength modes 
increases due to the presence of a large power reservoir 
in the form of the turbulent motions of the plasma in 
the same wavelength range. By contrast, in the case of 
helical fields, the power on large scales grows due to the 
transfer of power from the shorter wavelength modes. 

Apart from the transfer of power from the plasma mo- 
tions to the magnetic field at large scales, another effect 
of evolution of plasma and magnetic field perturbations 
is the turbulent decay of the power at short length scales 
due to the phenomenon of the direct turbulent cascade. 
At a given moment of time 77, this phenomenon leads 
to suppression of power in velocity and magnetic field 
modes on scales smaller than the size of the largest pro- 
cessed eddy. The size of the largest processed eddy is 
determined by the condition 



6f ~ Cm ~ v K (£ K )r), 



(8) 



where Vl(C) ls the characteristic velocity of the plasma 
motions on the scale £. From the definition of the kinetic 
energy power spectrum through the velocity two-point 



correlation Pk (k) 



\v K {L~l/k)\ 2 /k 2 one finds 



By a Kolmogorov-type spectrum we simply mean a fc — 5 / 3 spec- 
trum and ignore anisotropics that are known to exist in non- 
helical MHD [43j] . Such a spectrum can be derived from a phe- 
nomenological approach too; see 5.3.2. of Ref. [29ll , Our numeri- 
cal simulations 23 , 44] confirmed the Kolmogorov-type spectrum 
for a wide range of magnetic Prandtl numbers 



4 In fact, the two Reynolds numbers in the Universe are high 
enough to ensure the validity of the Kolmogorov-type phe- 
nomenological approach [30l ] ■ according to which there is self- 
similarity between kinetic and magnetic energy densities evolu- 
tion (see Sec. 7.3.4 of Ref. [H), i.e. 

£ M ~£ K ~£, (7) 
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that characteristic velocity vk on the distance scale L ~ 
fc -1 is v 2 K ~ k 2 Pj^. Since the power spectrum of plasma 
perturbations is Pk ~ re , we have ujf ~ fc ~ £ , so the 
above equation has the solution 



1/2 



(9) 



The energy of plasma perturbations on the scale £k is 
Ek = ^irk 2 P K ~ fc 2 ~ Cif 2 ~ Since the magnetic 

field on the scale ~ £a/ is in equipartition with the 
plasma, the energy density of the magnetic field evolves 
with time as 



E 



M 



E 



so that the strength of magnetic field evolves as 
#(eff) 



(10) 



(11) 



Below we demonstrate numerically (see Sec. Ill) that 
the evolution laws £m ~ rj L ' 2 ,EM ~ V 1 are indeed re- 
alized in the free turbulence decay regime. It should be 
noted that the "universality" of the decay law £j\/ oc r\ x l 2 
is not realized if we were to consider the magnetic field 
evolution separately from the velocity field evolution [2!| ; 
see also Refs. [HI, [H| for the magnetic field decay laws 
in the cosmological context. Accounting for the Loitsian- 
skii invariant for turbulence leads to the decay laws being 
dependent on the spectral shape; see Ref. [34| and ref- 
erences therein. It has also been claimed that the decay 
laws in the case of non-helical magnetic fields strongly de- 
pend on the initial conditions and can be different even 
when the helicity is extremely small [29j . 



2. Helical fields 

In the case of helical fields the evolution of Em and 
£m is determined directly by the condition of the conser- 
vation of magnetic helicity, J A ■ Bd 3 x. High Reynolds 
numbers allow us to follow the Kolmogorov-type phe- 
nomenological approach given above: see Sec. 4.2.3 of 
Ref. (29[. Accounting for magnetic helicity conserva- 
tion £m£,m = const 5 and combining Eqs. and £m ~ 
£ 3 / 2 /e, which follows from the dimensionless analysis 
(based on the Kolmogorov-type approach), we get |3(| 



d£ 



M 



drj 



£ 



-5/2 



M 



(12) 



which leads to the decay laws £m oc r\ 2 / 3 and £ 



2/3 



U OC 



Below we demonstrate numerically (see Sec. IIII B|) 
the appearance of the £m ~ t 2 ' 3 ,EM ~ t~ 2 ' 3 laws in 
the evolution of helical magnetic fields in the free decay 
regime. 



C. Simulations setup 

To model the evolution of the magnetic field and fluid 
perturbations we solve the compressible equations with 
the pressure given by p = pc 2 , where p is the gas density 
and c s = l/y/3 is the sound speed for an ultra-relativistic 
gas. Following our earlier work [35| . we solve the govern- 
ing equations for the logarithmic density In p, the velocity 
v, and the magnetic vector potential A, in the form 



Dlnp 
Dry 
Dt; 

dA 

drj 



-V-v, (13) 
JxB-c 2 Vlnp + / visc , (14) 
v x B + f M + XV 2 A, (15) 



where D/D77 = d/drj + v ■ V is the advective derivative, 
/vise = v (y 2 v + |VV ■ v + Gj is the viscous force in 
the compressible case with constant kinematic viscosity 
v and Gi = 2SyVj lnp as well as Sy = ^(Vij + vj i) — 
^SijVk : k being the trace-free rate of strain tensor. Fur- 
thermore, J = V x B/4tt is the current density. 

We use the Pencil Code [HJ with a resolution of 512 3 
meshpoints. An important difference of our simulations 
setup from those of previous studies is in the treatment 
of the backreaction of fluid perturbations onto the mag- 
netic field. Such treatment is important, especially on 
large length scales, because the spectrum of velocity per- 
turbations follows a white noise (Ek oc fc 2 ) spectrum 
to large scales [2(|. Large-scale fluid perturbations af- 
fect the magnetic field evolution at the largest scales and 
could lead to a transfer of power to the large-scale modes 
of the magnetic field, an effect similar to the inverse cas- 
cade developing even in the case of non-helical fields; see 
also Sec. [Till 

Prior to the simulation of the magnetic field decay we 
inject magnetic energy into the computational domain 
at scales corresponding to the phase transition eddy size 
see Ref. [3|| for more details. We approximate the mag- 
netic field injection by a delta function, allowing it to 
interact through MHD processes with a rest plasma. Af- 
ter several turnover times the initial sharp peak of the 
magnetic field starts to disappear and the magnetic field 
begins to spread over a wide range of the wavenumbers; 
see Fig. 4 of Ref. [35j]. In a few turnover times the spec- 
trum becomes established with a cut-off at small length 
scales and a well defined Kolmogorov-like integral scale 
fcjvf, and a fc 4 spectral shape at lar ge l ength scale. In 
contrast to the studies of Refs. [32l [33j we recover the 
spectral shape of the magnetic field at large scale. 6 This 
is in perfect agreement with previous analytical results 



Accounting for magnetic helicity conservation and assuming that 
the magnetic spectral energy is sharply peaked at £m a very 
rough estimate implies that i? eff oc Ca/^ 2 - 



6 Note that in Ref. [35j the spectral index between 3 and 4 due to 
a different choice of the system parameters and run time. 
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of Ref. [2l| based on causality and divergence free field 
arguments. 

We use simple power-law models for the decay of tur- 
bulence and scale the magnetic correlation length and 
magnetic field with temperature as follows: 



A = 



(16) 



(17) 



where B+ and T+ are the effective values of the magnetic 
field at the moment of generation and the temperature 
of the phase transition, respectively. Hence, the values of 
the parameters and tie describe the turbulent decay 
laws that differ from each other in the non-helical and 
helical cases. 

Qualitative arguments presented above show that the 
expected values for (n^, He) for non-helical and helical 
fields are (1/2,-1/2) and (2/3,-1/3), respectively. Be- 
low we show that this is indeed the case. 



III. SIMULATIONS RESULTS 

In Sec. II we have presented a phenomcnological de- 
scription of the scaling laws for the magnetic correlation 
length and the magnetic energy in the free turbulence 
decay regime. Below we address the same scaling laws 
based on our simulations. We also briefly review the re- 
sults of previous works. 
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FIG. 1: (color online). £m(£) for helical (thin, red) and non- 
helical (thick, blue) cases. 
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A. Non-helical Magnetic Field Evolution 



FIG. 2: (color online). £m(£) (solid) and £k(£) (dashed) for 
the helical (thin, red) and non-helical (thick, blue) cases. 



The scaling laws for the non-helical magnetic field evo 
lution have been studied throu gh d ifferent simulations by 



different groups; see Refs. 




46J and refer- 
the magnetic 



ences therein. As is stated in Ref. [2 
decay laws for the non-helical case strongly depend on 
initial conditions, and result in exponents n in the decay 
law £m(v) x V " that vary in the range 1.3 > n > 0.65. 
Note that the numerical and phenomenological studies 
performed in Refs. [HI, HH lead to £m(v) k V° A an d 
£m{v) V f° r a white noise spectrum, and this is 
in good agreement with the grid turbulence description 
of hydrodynamic turbulence [3?} ■ On the other hand, 
the 3D MHD simulations of Refs. [H] and the phe- 



nomenological study of Ref. [31[ show a slightly faster 
growth of the correlation length £m(?/) ^ V 1 ^ 2 with a 
magnetic energy decaying as £m °c t . The difference 
between two different scaling laws is probably due to dif- 
ferent initial conditions. In particular, the initial velocity 
field has traditionally been taken to be zero. By con- 
trast, here we have taken as initial condition the result 
of a self-consistent magnetically driven turbulence simu- 
lation. The numerical simulations of Ref. [441 show that 



the growth of the correlation length is almost indepen- 
dent of the magnetic Prandtl number with £m cx rj 1 ^ 2 
(Fig. [T]), while the exponent of the total magnetic en- 
ergy density decay is compatible with —1 (here closest 
to —0.9; Fig. [2]). Accounting for T cx a™ 1 and cx 
during the radiation dominated epoch, and the magnetic 
field strength = VS^Em we get the scaling indices 

for the decay of non- helical turbulence: = 1/2 and 
he = —1/2. As expected, for k <C fco early times, we 
find £ , M(fc,C) oc fc 4 and Eic(k,£) cx fc 2 ; see Figure [31 
Furthermore, even in this non-helical case the spectral 
energies increase with time for k <C kg, while for k 3> fco 
they decrease. 



B. Helical Magnetic Field Evolution 

As we have noted above, the presence of magnetic he- 
licity results in the development of an inverse cascade 
during which the correlation length is increasing while 
the total magnetic energy decreases. Similar to the non- 



7 




0.1 1.0 10.0 0.1 1.0 10.0 



FIG. 3: (color online). Evolution of Enj(k,(;) (solid) and 
E K {k,£) (dashed) versus k for 77 = 5, 10, 20, 50, and 100 
for the non-helical run. Thick lines are for 77 = 10. The red 
dash-dotted lines give the k 2 and k 4 scalings for comparison. 



FIG. 4: (color online). Evolution of E M (k,£,) (solid) and 
E K (k,£) (dashed) versus k for 77 = 5, 10, 20, 50, and 100 for 
the helical run. Thick lines are for rj = 10. The red dash- 
dotted lines give the k 2 and k 4 scalings for comparison. 



helical case there are basically two different approaches: 
(i) Refs. (3ol . l32l l33j assume exact conservation of mag- 
netic helicity and the magnetic field is the dominant con- 
tribution to the total energy density, i.e., £k/£m <C 1 
(where Ek is the total kinetic energy density of turbu- 
lence); (ii) other approaches are given in Refs. [27l. [3l|. 
In particular, Ref. [27] refers to a more general case with 
magnetic helicity evolving as H-uij]) oc rf 2s . Also, the 
ratio between kinetic and magnetic energy densities has 
in some studies assumed to be around 1; Ref. [3l[ assumes 
that the magnetic field evolves toward a force-free regime 
with constant magnetic helicity and with a constant ratio 
between magnetic and kinetic energy densities. All mod- 
els [28|, |30| - [33| show that the scaling laws are independent 
of the initial magnetic field spectrum. In fact, there are 
two main behaviors described: (i) Refs. [30l I32I |33[ claim 
£m(v) oc t?" 2 / 3 and £ M (rf) oc r? 2 / 3 ; (ii) RefsHil HJ 
claim ^m{v) 04 V 1 ^ 2 ? with £m(v) oc V 1 ^ 2 - I 11 both scal- 
ing laws £,m£m ~ const. The main difference between 
these two scaling laws consists in choosing the turbulence 
model. Refs. [13, Hl[ assume a force-free development 
of the MHD turbulence decay, while Ref. [H[ , see their 
Eq. (4), assumes a linear dependence between vorticity 
and Lorentz force. Our new numerical results support 
the former scenario (i). Fig. [3] shows the evolution of ki- 
netic and magnetic spectral energies. As we can see, the 
k 2 and fc 4 laws are established at large scales for kinetic 
and magnetic spectral energies, respectively. We can also 
see the slight increase of power at large scales, even in the 
case of non-helical fields. 

We have performed a study of the large-scale decay 
of a maximally helical magnetic field under conditions 
similar to those in the non- helical case, and for different 
magnetic Prandtl numbers as well as different values of 
magnetic resistivity. We have recovered the 2£m(&) oc fc 4 
spectral shape and the scaling laws as £m oc 77 2 / 3 and 
£m oc rj~ 2 / 3 , so that = 2/3, tie = —1/3 [47| : see also 



[29l | for more general discussion; see Figs. [T] and [2] Again 
these scaling laws are valid when the correlation length 
is greater than the damping scale, so that dissipation 
does not play an important role. Similarly to the non- 
helical run, we find for k -C fco and early times that 
Em{k,i) oc k A and E K {k,^) oc k 2 ; see Figure [4] In 
this case there is a strong inverse cascade with a strong 
increase of spectral energies with time for k <C kg. We 
can also see the constant magnetic power while the peak 
is moving toward large scales (inverse cascade). This 
corresponds to the constant helicity case. 



IV. IMPLICATION FOR COSMOLOGICAL 
EVOLUTION OF MAGNETIC FIELD 

Decay of cosmological MHD turbulence occurs to- 
gether with the cooling of the Universe and the increase of 
magnetic correlation length. Therewith, the correlation 
length increases only up to the point when the Universe 
reaches the temperature T = 1 eV [36[ . The initial values 
of correlation length and magnetic field strength, £0 an d 
Bq, at the temperature of magnctogenesis T*. together 
with the two scaling indices and tie, fully determine 
the large-scale magnetic field decay, and as a result the 
final configuration of the magnetic field. 

Phcnomcnological arguments as well as numerical sim- 
ulations show that ri£ = 1/2 and ue = —1/2 in non- 
helical case, while = 2/3 and tie = —1/3 in the helical 
case during the turbulent regime. The speed of growth of 
£m is constant, independently of the relation between £m 
and £d. This implies that, in the case in which the initial 
correlation length of magnetic field is comparable to the 
size of cosmological horizon at the epoch of magnetogen- 
esis, the final correlation length reaches 2 x 10~ 4 Mpc and 
6 x 10~ 3 Mpc for non-helical magnetic fields generated 
at EWPT and QCDPT, respectively. In the helical case, 
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FIG. 5: (color online). The correlation length £m (top 
panel) and the maximal allowed B max (bottom panel) for 
a primordial helical (thin red) and non-helical (thick blue) 
magnetic fields generated during QCDPT. Constraints on 
the magnetic field at T = 0.25eV are set to B max — 
4.5 x 10 _11 Gauss (£m = 2 x l(T 3 Mpc) and B max = 1.3 x 
l(T 9 Gauss (£ M = 5.6 x l(T 2 Mpc) for non-helical and helical 
cases, respectively. Dashed lines show areas where damping 
processes may reduce ideal estimates. 



FIG. 6: (color online). The correlation length £m (top panel) 
and the maximal allowed B max (bottom panel) for a primor- 
dial helical (thin red) and non-helical (thick blue) magnetic 
fields generated during EWPT. Constraints on the magnetic 



field at T = 0.25eV are set to B n 



= 1.3 x lO -1 " Gauss 
= 1.1 x l(T 10 Gauss 



(fM = 3.5 x l(T b Mpc) and B max 

(Cm = 3.1 x l(T 4 Mpc) for non-helical and helical cases, re- 
spectively. Dashed lines show areas where damping processes 
may reduce ideal estimates. 



the correlation length reaches 10 kpc for fields generated 
at the QCDPT and 0.1 Mpc for the EWPT. 

The evolutionary paths of magnetic field strength and 
correlation length in time and in the -B,£m parameter 
space are shown in Figs. [5H7J In principle, the free tur- 
bulence decay periods in the early Universe are intermit- 
tent, with periods of viscously damped evolution (32[. So, 
starting from 10 MeV, Figs. [5] and |5] show only maximal 
values of helical and non-helical magnetic fields and their 
correlation length from QCD and EW phase transitions 
without accounting for any possible damping processes 
that can affect the given scaling laws. These maximal 
values are derived using several assumptions: the mag- 
netic correlation length during phase transition matches 
the bubble size, and magnetic fields are excited with max- 
imal amplitudes allowed by the BBN limit. In fact, be- 
fore neutrino decoupling the viscous damping force f v i sc 
in Eq. (14) grows as i^Z m fp, v ~ r/ 4 , where Z m f Pj „ is the 
neutrino mean free path. This growth is much faster 
than the £m ~ rj 1 ^ 2 growth of the integral scale of the 
magnetic field. This means that, even if £m 3> ^ m fp, v 
at the moment of magnetogenesis, l m f p catches up with 
£a/ at a later time ?y V i SC . Starting from this time and 
up to the moment of neutrino decoupling, the magnetic 
field stops to decay, B ~ const, because the fluid mo- 
tions are damped by viscosity, v ~ 0, so that there is 
no coupling of magnetic field to the fluid in this regime. 



However, turbulence re-starts after the neutrino decou- 
pling, so that the system returns to the same evolution- 
ary track shown in Fig. [7] just after neutrino decoupling. 
At lower temperatures, when the viscosity is provided by 
photon streaming, the viscous damping scale grows as 
v ~ ^mfp. 7 ~ "H 2 , where / mtp . 7 is the photon mean free 
path. In this time interval the growth of v is again faster 
than the growth of £m , so that the episode of viscously 
damped evolution repeats when £ m f Pj 7 reaches £m- This 
could again delay the advance of (B,^b) along the evo- 
lutionary track shown in Fig. UJ The end point of the 
evolutionary track at the end of the radiation-dominated 
era is well defined by the condition that the correlation 
length of magnetic field should not be shorter than the 
Silk damping scale times the Alfven velocity 32(. The 
loci of the possible end points of the evolution are shown 
by the inclined thick solid (green/orange) lines in Fig.[7J 



V. CONCLUSION 

Our study shows that magnetic fields generated during 
phase transitions are comparable with the observational 
lower bound even if we account for large-scale decay as 
well as additional Alfven wave-induced damping. The 
extremely low values of the smoothed magnetic field [25| 
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FIG. 7: (color online). Cosmological evolution of B^ c ' and 
£m for magnetic fields generated at the EWPT (green) and 
QCDPT (orange). Arrows show the evolution of the strength 
and integral scale of helical and non-helical fields during the 
radiation dominated era up to their final values. Thick solid 
line(s) show possible present day strength and integral scale 
of the phase-transition generated magnetic fields. 

do not imply that the effective magnetic field in the range 
I pc-1 kpc are small enough to result in observational 



changes in blazar emission spectra. The advantage of 
using the effective magnetic field lies in its independence 
of the spectral shape. In summary, if the magnetic field 
has been generated during a phase transition, its corre- 
lation length is strongly limited. If future observations 
were to detect a weak magnetic field < 10~ 14 -10~ 15 G 
with a typical correlation length of the order of a few 
pc, this could serve as an indication of magnctogenesis 
during EWPT, while a somewhat stronger field with a 
correlation length of the order of kpc might indicate the 
presence of QCDPT magnctogenesis. 
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